Skin transcriptome analysis in Brangus cattle under heat stress.

Heat stress is a major factor that negatively affects animal welfare and production systems. Livestock should adapt to tropical and subtropical areas and to meet this, composite breeds have been developed. This work aimed to evaluate gene expression profiles in the skin of Brangus cattle under heat stress using a case-control design, and to correlate this with skin histological characteristics. Two groups of bulls were set using rectal temperature as a criterion to define stress conditions: stressed (N = 5) and non-stressed (N = 5) groups. Skin transcriptomics was performed and correlations between breed composition, phenotypic and skin histological traits were evaluated. Results showed 4309 differentially expressed genes (P < 0.01), 2113 downregulated and 2196 upregulated. Enrichment and ontology analyses revealed 132 GO terms and 67 pathways (P < 0.01), including thermogenesis, glycolysis, gluconeogenesis, mitochondrial activity, antioxidant and immune response, and apoptosis. The identity of the terms and pathways indicated the diversity of mechanisms directed to relieve the animals' suffering, acting from simple passive mechanisms (conduction, convection and radiation) to more complex active ones (behavioural changes, evaporation, vasodilation and wheezing). Furthermore, significant differences between phenotypic and skin histological traits and correlations between pairs of traits suggested a direction towards heat dissipation processes. In this sense, number of vessels was positively correlated with number of sweat glands (P < 0.001) and both were positively correlated with zebuine genetic content (P < 0.05 and P < 0.01, respectively), gland size was positively correlated with epidermal thickness and negatively with hair length (P < 0.05), and epidermal thickness was negatively correlated with gland-epidermis distance (P < 0.0005). These results support the notion that response to heat stress is physiologically complex, producing significant changes in the expression of genes involved in several biological pathways, while the animal's ability to face it depends greatly on their skin features.

Exploring genetic diversity and variation of Ovar-DRB1 gene in Sudan Desert Sheep using targeted next-generation sequencing.

INTRODUCTION: The Ovar-DRB1 gene, a crucial element of the Major Histocompatibility Complex (MHC) Class II region, initiates adaptive immunity by presenting antigens to T-cells. Genetic diversity in sheep, particularly in MHC Class II genes like Ovar-DRB1, directly influences the specturm of presented antigens impacting immune responses and disease susceptability. Understanding the allelic diversity of Ovar-DRB1 gene in Sudan Desert Sheep (SDS) is essential for uncovering the genetic basis of immune responses and disease resistance, given the the breeds significance in Sudan's unique environment. METHODS: Utilizing Targeted Next-Generation Sequencing (NGS) we explore allelic diversity in Ovar-DRB1 gene within SDS. Successfully ampliying and and sequencing the second exon of this gene in 288 SDS samples representing six breeds provided a comprehensive allelic profile, enabling a detalied examination of the gene's genetic makeup. RESULTS: We identifed forty-six alleles, including four previously unreported, enrichness the genetic diversity of SDS breeds. These alleles exhibiting non-uniform distribution, varying frequencies across breeds, indicating a breed-specific genetic landscape. Certain alleles, known and novel, show higher frequencies in specific populations, suggesting potential associations with adaptive immune responses. Identifying these alleles sets the stage for investigating their functional roles and implications for disease resistance. Genetic differentiation among SDS breeds, as indicated by FST values and clustering analyses, highlights a unique genetic makeup shaped by geographic and historical factors. These differentiation patterns among SDS breeds have broader implications for breed conservation and targeted breeding to enhance disease resistance in specific populations. CONCLUSION: This study unveils Ovar-DRB1 gene allelic diversity in SDS breeds through targeted NGS and genetic analyses, revealing new alleles that underscore the breeds' unique genetic profile. Insights into the genetic factors governing immune responses and disease resistance emerge, promising for optimization of breeding strategies for enhanced livestock health in Sudan's unique environment.

Identification of Anaplasma marginale, Babesia bovis and Babesia bigemina resistance alleles in Crioulo Lageano cattle using PCR-SBT and BoLA-DRB3 gene sequencing.

Introduction: The BoLA-DRB3 gene in cattle is associated with tolerance to several infectious diseases, such as neosporosis, dermatophilosis, leukosis, and mastitis. Methods: This study used PCR-SBT and BoLA-DRB3 gene sequencing to determine the association between the presence or absence of Anaplasma marginale, Babesia bovis, and Babesia bigemina infections in 208 Crioulo Lageano cattle and alleles present in the population. The chi-square test and odds ratio analysis were employed to establish the association. Results: Of the BoLA-DRB3 gene alleles present in the population, two alleles were significantly associated with resistance to A. marginale infections: BoLA-DRB3001:01 (p < 0.001; OR = 0.224), which had a frequency of 7.93%, and BoLA-DRB3024:06 (p = 0.007; OR < 0.00001), which had a frequency of 0.72%. Regarding B. bovis infection, the BoLA-DRB3*011:01 allele (p = 0.002; OR = 0.271) had a frequency of 6% in the population and was associated with resistance to the infection. None of the alleles was associated with resistance to infection by B. bigemina. Discussion: The Crioulo Lageano breed has alleles that may confer resistance against infection by A. marginale and B. bovis.

Phenotypic evaluation of genetic resistance to the tick Rhipicephalus (Boophilus) microplus in Argentine Creole cattle.

The objective of this work was to characterize the Argentine Creole cattle breed through the identification of individual phenotypic variations in the levels of infestation with Rhipicephalus (Boophilus) microplus. We evaluated 179 heifers exposed to successive artificial infestations from 2015 to 2018, achieving a total of 663 observations. Tick counts were assessed with the linear mixed model, considering year of evaluation, time of infestation, dam's age and nutritional status during the evaluated period as fixed effects. The average tick count value obtained allowed to classify the breed as highly resistant to the tick charge (99.3%). Although the previous nutritional condition of the animals did not affect the individual charge response, weight gain during the trial showed a significantly negative correlation. We conclude that the Argentine Creole breed is an attractive genetic alternative for cattle breeding in endemic regions, either as a pure breed or a cross-breed.

Preliminary transcriptomic analysis of peripheral blood from German Shepherd dogs with degenerative joint disease for the identification of diagnostic biomarkers.

OBJECTIVE: Degenerative joint disease (DJD) includes a group of disorders characterised by the deterioration of the articular cartilage. In this study, we investigated the transcriptomic profile of peripheral blood in German Shepherd dogs with DJD to identify putative diagnostic biomarkers. METHODS: Differential gene expression (DGE) and gene ontology (GO) analyses of the bulk RNA-seq experiment were performed in a cohort of 12 adult dogs (five cases and seven controls, classified by clinical and radiographic analyses). RESULTS: Radiographs of cases revealed severe signs of progressive DJD. Two up-regulated (LOC106559672 and THBS4) and one down-regulated (LOC106559235) differentially expressed genes (adjusted p value < 0.05) were identified. The DGE with log2 fold change < -1.5 and > 1.5 and non-adjusted p < 0.01 were selected for GO analysis. No significant enrichment terms were observed in the selected threshold. CONCLUSION: The gene-encoding protein THBS4 is correlated with DJD severity and long noncoding RNA LOC106559235 is probably involved in the DJD process. The THBS4 gene should be considered a good biomarker for DJD in dogs. Future studies using independent cohorts will be necessary to validate the present results.

Coat score. A possible explanation for the zebuine selective sweep located on bovine chromosome 5: 47,670,001-48,100,000 bp.

Over 65% of the world's cattle population resides in warm areas where heat stress conditions limit the breed of European taurine cattle. Composite breeds were developed to retain the main traits of both parental breeds. The skin plays a central role in animal response to heat stress. Research on the genetic architecture of skin traits has identified genes and regions related to warm resistance skin features. The aim of this study was to determine whether the indicine proportion accounted for coat type or whether there were genes of large effect segregating in Brangus. Bulls (n = 108) were genotyped using microarrays and their coat score and hair length were evaluated. Indicine-taurine genome-wide composition was estimated and GWAS was performed. Although significant correlations between indicine proportion and traits were not observed, four windows of SNPs on BTA4 and BTA5 explained more than 2% of the trait variance. The GWAS for coat score in summer showed the main peak on BTA5:46,941,446-48,030,219 bp, accounting for 4.65% of the variance. Our results suggest that the variation in coat score and undercoat hair length in Argentinian Brangus bulls is associated with the presence of some particular gene variants, rather than with the whole indicine genetic content.

Characterization of experimental cattle populations from Argentina with a low-density SNP genotyping panel.

The objectives of the present experiment were to evaluate a low-density SNP array designed for the molecular characterisation of gene banks and to assess the genetic diversity and population structure of beef cattle herds from an Argentinean research station. Forty-nine animals from three breeds (Angus, Hereford, and Argentinean Creole) were genotyped using the multi-species IMAGE001 60-K SNP array (10 K for cattle). Genotypes of other 19 cattle populations from Argentina, other American countries, and Europe were included in the study. Of special interest was the characterization of the Argentinean Creole, the only autochthonous cattle breed in the country. Due to the merging of different datasets, approximately 5 K SNPs were effectively used. Genetic differentiation (FST ), principal component analysis, neighbour-joining tree of Reynolds distances and ancestry analysis showed that autochthonous American breeds are clearly differentiated, but all have genetic influences of Iberian cattle. The analysed herds of Argentinean Creole showed no evidence of recent admixture and represent a unique genetic pool within local American breeds. An experimental herd and the local commercial Hereford population have also diverged, probably due to the influence of current selection objectives in the breed. Our results illustrate the utility of using low-cost, low density SNP arrays in the evaluation of animal genetic resources. This type of panels could become a very useful resource in developing countries, where most endangered cattle breeds are located. The results also reinforce the importance of experimental herds as reservoir of genetic diversity, particularly in the case of local breeds under-represented in traditional production systems.

Genome-wide scan for signatures of selection in the Brangus cattle genome.

Brangus is a composite cattle breed developed with the objective of combining the advantages of Angus and Zebuine breeds (Brahman, mainly) in tropical climates. The aim of this work was to estimate breed composition both genome-wide and locally, at the chromosome level, and to uncover genomic regions evidencing positive selection in the Argentinean Brangus population/nucleus. To do so, we analysed marker data from 478 animals, including Brangus, Angus and Brahman. Average breed composition was 35.0% ± 9.6% of Brahman, lower than expected according to the theoretical fractions deduced by the usual cross-breeding practice in this breed. Local ancestry analysis evidenced that breed composition varies between chromosomes, ranging from 19.6% for BTA26 to 56.1% for BTA5. Using approaches based on allelic frequencies and linkage disequilibrium, genomic regions with putative selection signatures were identified in several chromosomes (BTA1, BTA5, BTA6 and BTA14). These regions harbour genes involved in horn development, growth, lipid metabolism, reproduction and immune response. We argue that the overlapping of a chromosome segment originated in one of the parental breeds and over-represented in the sample with the location of a signature of selection constitutes evidence of a selection process that has occurred in the breed since its take off in the 1950s. In this regard, our results could contribute to the understanding of the genetic mechanisms involved in cross-bred cattle adaptation and productivity in tropical environments.

BoLA-DRB3 gene haplotypes show divergence in native Sudanese cattle from taurine and indicine breeds.

Autochthonous Sudanese cattle breeds, namely Baggara for beef and Butana and Kenana for dairy, are characterized by their adaptive characteristics and high performance in hot and dry agro-ecosystems. They are thus used largely by nomadic and semi-nomadic pastoralists. We analyzed the diversity and genetic structure of the BoLA-DRB3 gene, a genetic locus linked to the immune response, for the indigenous cattle of Sudan and in the context of the global cattle repository. Blood samples (n = 225) were taken from three indigenous breeds (Baggara; n = 113, Butana; n = 60 and Kenana; n = 52) distributed across six regions of Sudan. Nucleotide sequences were genotyped using the sequence-based typing method. We describe 53 alleles, including seven novel alleles. Principal component analysis (PCA) of the protein pockets implicated in the antigen-binding function of the MHC complex revealed that pockets 4 and 9 (respectively) differentiate Kenana-Baggara and Kenana-Butana breeds from other breeds. Venn analysis of Sudanese, Southeast Asian, European and American cattle breeds with 115 alleles showed 14 were unique to Sudanese breeds. Gene frequency distributions of Baggara cattle showed an even distribution suggesting balancing selection, while the selection index (ω) revealed the presence of diversifying selection in several amino acid sites along the BoLA-DRB3 exon 2 of these native breeds. The results of several PCA were in agreement with clustering patterns observed on the neighbor joining (NJ) trees. These results provide insight into their high survival rate for different tropical diseases and their reproductive capacity in Sudan's harsh environment.

Epidemiology of Cranial Cruciate Ligament Rupture and Patellar Luxation in Dogs from the Province of Buenos Aires, Argentina.

OBJECTIVE: The aim of this study was to investigate the prevalence of cranial cruciate ligament rupture and patellar luxation and the associated risk factors in dogs. MATERIALS AND METHODS: A total of 13,072 clinical records of dogs were reviewed from School Hospital (Faculty of Veterinary Sciences, National University of La Plata). Data of age, breed, sex, body weight, patellar luxation and cranial cruciate ligament rupture condition were registered. Chi-squared and Fisher's exact tests were used to compare the prevalence of cranial cruciate ligament rupture and patellar luxation with the variables and then univariable logistic regression was used to evaluate the risk of having cranial cruciate ligament rupture and patellar luxation. Multivariable logistic regression was used including all variables to assess the odds of having patellar luxation and cranial cruciate ligament rupture. RESULTS: Of 13,072 patients treated, 72 and 51 had cranial cruciate ligament rupture and patellar luxation respectively. Sex was not a major risk factor for either condition. Adult (odds ratio [OR] = 8.2) and senior (OR = 4.3) patients had increased risk of having cranial cruciate ligament rupture, while for patellar luxation age was not a risk factor. Groups 2, 3 and 8 were more likely to have cranial cruciate ligament rupture (OR = 5.5, OR = 9.1 and OR = 2.6), and group 11 had lower risk of having patellar luxation (OR = 0.08). Maxi (OR = 2.4) and giant (OR = 6.0) breeds had higher risk of having cranial cruciate ligament rupture, and medium and maxi breeds had higher risk of patellar luxation (OR = 0.05 and OR = 0.3). Multivariate OR test confirmed that age (adult), body size (giant and maxi) and breed group (Group 3) were significantly associated with having cranial cruciate ligament rupture, and age was associated with having patellar luxation. CLINICAL SIGNIFICANCE: This is the first epidemiological study of cranial cruciate ligament rupture and patellar luxation in dogs from School Hospital (Faculty of Veterinary Sciences, National University of La Plata). Giant and large adult dogs from the Molossoid and Terrier breeds were more likely to have cranial cruciate ligament rupture, while mixed and large dog breeds showed the lowest risk of having patellar luxation.

BoLA-DRB3 genetic diversity in Highland Creole cattle from Bolivia.

The genetic diversity of the BoLA-DRB3 gene has been reported in different cattle breeds owing to its central role in the immune response. However, it is still unknown in hundreds of cattle breeds, especially native populations. Here, we studied BoLA-DRB3 genetic diversity in Highland Creole cattle (CrAl) from Western Bolivia, raised at altitudes between 3800 and 4200 m. DNAs from 48 CrAl cattle were genotyped for BoLA-DRB3 exon 2 alleles using polymerase chain reaction-sequence-based typing (PCR-SBT). The results were compared with 1341 previously reported data from Tropical Creole cattle and other breeds raised in the region. Twenty-three BoLA-DRB3 alleles were identified in CrAl, including the BoLA-DRB3*029:02 variant previously detected in other Creole cattle. Observed and expected heterozygosity were 0.87 and 0.93, respectively. Nucleotide diversity and the number of pairwise difference values were 0.078 and 19.46, respectively. The average number of nonsynonymous and synonymous substitutions were 0.037 and 0.097 for the entire BoLA-DRB3 exon 2, and 0.129 and 0.388 for the antigen-binding site, respectively. Venn analysis and the review of the IPD-MHC database and the literature showed that 2 of 64 alleles were only detected in CrAl, including BoLA-DRB3*029:01 previously reported in African cattle and *048:01 detected in Philippine cattle. Two additional alleles, BoLA-DRB3*007:02 and *029:02, were only present in CrAl and Lowland Creole cattle. Principal Component Analysis (PCA) showed that Bolivian Creole cattle breeds were closely located but they were distant from the Colombian Hartón del Valle Creole. FST analysis showed a low degree of genetic differentiation between Highland and Lowland Bolivian Creole cattle (FST = 0.015). The present results contribute to increasing our knowledge of BoLA-DRB3 genetic diversity in cattle breeds.

Author Correction: The genetic ancestry of American Creole cattle inferred from uniparental and autosomal genetic markers.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Characterization of bovine MHC DRB3 diversity in global cattle breeds, with a focus on cattle in Myanmar.

BACKGROUND: Myanmar cattle populations predominantly consist of native cattle breeds (Pyer Sein and Shwe), characterized by their geographical location and coat color, and the Holstein-Friesian crossbreed, which is highly adapted to the harsh tropical climates of this region. Here, we analyzed the diversity and genetic structure of the BoLA-DRB3 gene, a genetic locus that has been linked to the immune response, in Myanmar cattle populations. METHODS: Blood samples (n = 294) were taken from two native breeds (Pyer Sein, n = 163 and Shwe Ni, n = 69) and a cattle crossbreed (Holstein-Friesian, n = 62) distributed across six regions of Myanmar (Bago, n = 38; Sagaing, n = 77; Mandalay, n = 46; Magway, n = 46; Kayin, n = 43; Yangon, n = 44). In addition, a database that included 2428 BoLA-DRB3 genotypes from European (Angus, Hereford, Holstein, Shorthorn, Overo Negro, Overo Colorado, and Jersey), Zebuine (Nellore, Brahman and Gir), Asian Native from Japan and Philippine and Latin-American Creole breeds was also included. Furthermore, the information from the IPD-MHC database was also used in the present analysis. DNA was genotyped using the sequence-based typing method. DNA electropherograms were analyzed using the Assign 400ATF software. RESULTS: We detected 71 distinct alleles, including three new variants for the BoLA-DRB3 gene. Venn analysis showed that 11 of these alleles were only detected in Myanmar native breeds and 26 were only shared with Asian native and/or Zebu groups. The number of alleles ranged from 33 in Holstein-Friesians to 58 in Pyer Seins, and the observed versus unbiased expected heterozygosity were higher than 0.84 in all the three the populations analyzed. The FST analysis showed a low level of genetic differentiation between the two Myanmar native breeds (FST = 0.003), and between these native breeds and the Holstein-Friesians (FST <  0.021). The average FST value for all the Myanmar Holstein-Friesian crossbred and Myanmar native populations was 0.0136 and 0.0121, respectively. Principal component analysis (PCA) and tree analysis showed that Myanmar native populations grouped in a narrow cluster that diverged clearly from the Holstein-Friesian populations. Furthermore, the BoLA-DRB3 allele frequencies suggested that while some Myanmar native populations from Bago, Mandalay and Yangon regions were more closely related to Zebu breeds (Gir and Brahman), populations from Kayin, Magway and Sagaing regions were more related to the Philippines native breeds. On the contrary, PCA showed that the Holstein-Friesian populations demonstrated a high degree of dispersion, which is likely the result of the different degrees of native admixture in these populations. CONCLUSION: This study is the first to report the genetic diversity of the BoLA-DRB3 gene in two native breeds and one exotic cattle crossbreed from Myanmar. The results obtained contribute to our understanding of the genetic diversity and distribution of BoLA-DRB3 gene alleles in Myanmar, and increases our knowledge of the worldwide variability of cattle BoLA-DRB3 genes, an important locus for immune response and protection against pathogens.

Naturally monozygotic quadruplets in a Braford cow confirmed by DNA analysis: A case report.

Cattle are a monotocous species, despite naturally conceived multiple births are sometimes observed. Although the number of twins has consistently increased, triplet and quadruplet pregnancies represent 0.015% and 0.004% of the total births, respectively. Multiple births are the result of multiple ovulation and/or the spontaneous cleavage of one fertilized oocyte, which is known as monozygotic (MZ) twinning. In cattle, approximately 5% to 14% of all twin births are MZ, and births with more than two MZ calves are extremely rare. Monozygotic animals are genetically identical, and those derived from two or more zygotes are genetically different. Furthermore, the presence of placental vascular anastomosis can result in foetal chimerism. Notwithstanding, animals born as single calves can be chimeras when one of the foetal twins dies undetected in utero. Here, we used DNA testing to study the zygotic condition of an unusual female quadruplet born from a Braford cow bred in a multi-sire natural mating system without hormone stimulation. Two tissues with different embryological origin were sampled to test zygosity and possible chimerism. The results showed an identical genotype, confirming they all originated in an MZ pregnancy and suggesting the lack of chimerism in all animals. The use of MZ twins in breeding and selection systems provides an alternative to the conventional progeny testing. Some works have suggested a genetic background of MZ twins in humans. This female and her daughters could be the founders of a lineage to study the possible inheritance of MZ multiple births in cattle.

Application of longitudinal data analysis allows to detect differences in pre-breeding growing curves of 24-month calving Angus heifers under two pasture-based systems with differential puberty onset.

BACKGROUND: Longitudinal data analysis contributes to detect differences in the growing curve by exploiting all the information involved in repeated measurements, allowing to distinguish changes over time within individuals, from differences in the baseline levels among groups. In this research, longitudinal and cross-sectional analysis were compared to evaluate differences in growth in Angus heifers under two different grazing conditions, ad libitum (AG) and controlled (CG) to gain 0.5 kg day-1 . RESULTS: Longitudinal mixed models show differences in growing curve parameters between grazing conditions, that were not detected by cross-sectional analysis. Differences (P < 0.05) in first derivative of growth curves (daily gain) until 289 days were observed between treatments, AG being higher than CG. Correspondingly, pubertal heifer proportion was also higher in AG at the end of rearing (AG, 0.94; CG, 0.67). CONCLUSION: In longitudinal studies, the power to detect differences between groups increases by exploiting the whole information of repeated measures, modelling the relation between measurements performed on the same individual. Under a proper analysis, valid conclusion can be drawn with fewer animals in the trial, improving animal welfare and reducing investigation costs. © 2019 Society of Chemical Industry.

The genetic ancestry of American Creole cattle inferred from uniparental and autosomal genetic markers.

Cattle imported from the Iberian Peninsula spread throughout America in the early years of discovery and colonization to originate Creole breeds, which adapted to a wide diversity of environments and later received influences from other origins, including zebu cattle in more recent years. We analyzed uniparental genetic markers and autosomal microsatellites in DNA samples from 114 cattle breeds distributed worldwide, including 40 Creole breeds representing the whole American continent, and samples from the Iberian Peninsula, British islands, Continental Europe, Africa and American zebu. We show that Creole breeds differ considerably from each other, and most have their own identity or group with others from neighboring regions. Results with mtDNA indicate that T1c-lineages are rare in Iberia but common in Africa and are well represented in Creoles from Brazil and Colombia, lending support to a direct African influence on Creoles. This is reinforced by the sharing of a unique Y-haplotype between cattle from Mozambique and Creoles from Argentina. Autosomal microsatellites indicate that Creoles occupy an intermediate position between African and European breeds, and some Creoles show a clear Iberian signature. Our results confirm the mixed ancestry of American Creole cattle and the role that African cattle have played in their development.

Genomic structural diversity in Criollo Argentino horses: Analysis of copy number variations.

Copy number variation (CNV) has been proved to be widespread in human, animal and plant genomes. Together with single nucleotide polymorphisms (SNPs), CVNs play a key role in genetic diversity. In this study, genome-wide detection of CNVs was performed based on SNP data from 24 Criollo Argentino horses genotyped with the GGP Equine70k array. Overall, 165 CNVs meeting stringent quality control criteria were identified and then aggregated into 87 CNV regions (CNVRs), representing a horse genome coverage of 13.69 Mb. Functional analysis of CNVRs allowed the identification of 337 genes implicated in a wide range of biological functions such as signal transducer activity (olfactory receptors), receptor activities and binding. Furthermore, enrichment analysis showed that the most represented protein classes (over 25%) were immunoglobulin receptor subfamily, immunoglobulins and major histocompatibility complex antigen (beta-2-microglobulin). To the best of our knowledge, this is the first report of CNV in Criollo Argentino horses.

Runs of homozygosity in a selected cattle population with extremely inbred bulls: Descriptive and functional analyses revealed highly variable patterns.

The analysis of runs of homozygosity (ROH), using high throughput genomic data, has become a valuable and frequently used methodology to characterize the genomic and inbreeding variation of livestock and wildlife animal populations. However, this methodology has been scarcely used in highly inbred domestic animals. Here, we analyzed and characterized the occurrence of ROH fragments in highly inbred (HI; average pedigree-based inbreeding coefficient FPED = 0.164; 0.103 to 0.306) and outbred Retinta bulls (LI; average FPED = 0.008; 0 to 0.025). We studied the length of the fragments, their abundance, and genome distribution using high-density microarray data. The number of ROH was significantly higher in the HI group, especially for long fragments (>8Mb). In the LI group, the number of ROH continuously decreased with fragment length. Genome-wide distribution of ROH was highly variable between samples. Some chromosomes presented a larger number of fragments (BTA1, BTA19, BTA29), others had longer fragments (BTA4, BTA12, BTA17), while other ones showed an increased ROH accumulation over specific loci (BTA2, BTA7, BTA23, BTA29). Similar differences were observed in the analysis of 12 individuals produced by a similar inbred event (FPED3 = 0.125). The correlation between the fraction of the genome covered by ROH (FROH) and FPED was high (0.79), suggesting that ROH-based estimations are indicative of inbreeding levels. On the other hand, the correlation between FPED and the microsatellite-based inbreeding coefficient (FMIC) was only moderate (r = 0.44), suggesting that STR-based inbreeding estimations should be avoided. Similarly, we found a very low correlation (r = -0.0132) between recombination rate and ROH abundance across the genome. Finally, we performed functional annotation analyses of genome regions with significantly enriched ROH abundance. Results revealed gene clusters related to pregnancy-associated proteins and immune reaction. The same analysis performed for regions enriched with recently formed ROH (> 8 Mb) showed gene clusters related to flagellum assembly. In both cases, the processes were related to male and female reproductive functions, which may partially explain the reduced fertility associated with inbred populations.

Genetic diversity of BoLA-DRB3 in South American Zebu cattle populations.

BACKGROUND: Bovine leukocyte antigens (BoLAs) are used extensively as markers of disease and immunological traits in cattle. However, until now, characterization of BoLA gene polymorphisms in Zebu breeds using high resolution typing methods has been poor. Here, we used a polymerase chain reaction sequence-based typing (PCR-SBT) method to sequence exon 2 of the BoLA class II DRB3 gene from 421 cattle (116 Bolivian Nellore, 110 Bolivian Gir, and 195 Peruvian Nellore-Brahman). Data from 1416 Taurine and Zebu samples were also included in the analysis. RESULTS: We identified 46 previously reported alleles and no novel variants. Of note, 1/3 of the alleles were detected only in Zebu cattle. Comparison of the degree of genetic variability at the population and sequence levels with genetic distance in the three above mentioned breeds and nine previously reported breeds revealed that Zebu breeds had a gene diversity score higher than 0.86, a nucleotide diversity score higher than 0.06, and a mean number of pairwise differences greater than 16, being similar to those estimated for other cattle breeds. A neutrality test revealed that only Nellore-Brahman cattle showed the even gene frequency distribution expected under a balanced selection scenario. The FST index and the exact G test showed significant differences across all cattle populations (FST = 0.057; p <  0.001). Neighbor-joining trees and principal component analysis identified two major clusters: one comprising mainly European Taurine breeds and a second comprising Zebu breeds. This is consistent with the historical and geographical origin of these breeds. Some of these differences may be explained by variation of amino acid motifs at antigen-binding sites. CONCLUSIONS: The results presented herein show that the historical divergence between Taurine and Zebu cattle breeds is a result of origin, selection, and adaptation events, which would explain the observed differences in BoLA-DRB3 gene diversity between the two major bovine types. This allelic information will be important for investigating the relationship between the major histocompatibility complex and disease, and contribute to an ongoing effort to catalog bovine MHC allele frequencies according to breed and location.